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Equalious are derived from which the temperature dependence of both the specific heat and the thermal 
diffusivity of a spherical sample of material can be calculated from observations of the time dependence of 
the surface temperature and the time-rate of energy loss from the sample as it cools. The derivation takes 
into account the nonuniformity of the interior temperature field of the sample, and the reuniting equations 
can be applied not only to radiative cooling, but also to any other cooling mechanism that does not violate 
the assumed spherical symmetry. The analysis excludes change of phase, but it does take thermal expan- 
sion into account. To permit the making of estimates neeessary for the design of radiative cooling ex- 
periments, a universal temperature-time cooling curve is derived for the post-transient cooling regime of a 
radiating sphere of any size with arbitrary, but constant, thermal parameters. 
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1. Introduction 

The analysis presented in this paper is an outgrowth of 
a proposal made by J. H. Colwell [1,2] to determine the 
high-temperature values of the specific heat, thermal dif- 
fusivity, and total hemispherical emissivity of a spherical 
sample of refractory material by making independent 
optical observations of the surface temperature of the 
sample and its time-rate of energy loss as it cools by free 
radiation into a cold vacuum. The original proposal was 
made in the context of an experiment to be conducted on 
board the space shuttle, and envisaged induction heating 
of the sample. With this mode of heating, the total heat 
content of the spherical sample and its interior 
temperature field at the start of the observational run 
would be unknown. However, after an interval on the 
order of the characteristic thermal decay time of the sam- 
ple, the interior temperature field would settle into the 
"post-transient regime" in which the interior field would 
be entirely determined by the time-dependence of the 
surface-temperature. Thus, in the post-transient regime 
it should, in principle, be possible to determine the 
temperature dependence of the thermal parameters from 



*Centerfor Chemical Engineering, National Engineering Laboratory 



a knowledge of the time dependenee of the surface 
temperature and the rate of energy loss. The determina- 
tion of the emissivity is a trivial matter, since it is propor- 
tional to the rate of energy loss divided by the fourth 
power of the surface temperature- 
Stated in mathematical terms, the cooling sample 
could be regarded as mapping the temperature- 
dependent specific heat and diffusivity over into the 
time-dependent surface temperature and energy loss 
functions, and the analytical problem then consists of in- 
verting this mapping so as to be able to express the two 
unknown thermal parameters in terms of the two ob- 
served time functions. The first step in carrying out this 
inversion, the results of which are summarized in section 
2, is to find the "surface-driven solution" of the Fourier 
equation for specified temperature-dependent specific 
heat and thermal diffusivity. This solution is completely 
specified by the time dependence of the surface 
temperature, and the time rate of change of the total heat 
content can be calculated from it. If this calculated rate 
of change is then equated to the fourth power of the sur- 
face temperature in aeeordance with the Stefan- 
Boltzmann radiation law, a nonlinear ordinary differen- 
tial equation (of infinite order) results whieh ean be 
iteratively solved (in truncated form) for the ease of con- 
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stant thermal parameters to yield a universal 
temperature versus time dependence for the post- 
transient regime. This solution, which is presented in 
section 3, is useful for making the various estimates that 
are necessary for the design of a radiative-cooling 
calorimetry experiment. At the end of section 3 an in- 
tegral equation is given that could also be used as the 
basis of an iterative solution of the post-transient predic- 
tive problem. 

In section 4 the surface-driven solution that is sum- 
marized in section 2 is inverted so as to yield expressions 
for specific heat and diffusivity in terms of the observed 
time-dependent surface temperature and time-rate of 
energy loss. These expressions constitute the desired 
solution of the calorimetric problem. The thermal 
parameters are expressed both in terms of truncated ex- 
pansions whose coefficients involve higher-order time 
derivatives of the observed functions, as well as in terms 
of integral expressions involving retrospective weighted 
averages of the observed time-dependent functions. The 
truncated expansions, which are easier to apply than the 
integral expressions, ought to suffice for analyzing most 
post-transient experiments. In fact, in many experiments 
the simple approximate expressions given in eqs (44) and 
(45) will be sufficiently accurate. In section 4 an estimate 
of the range of validity of these simple expressions, as 
well as the range of validity of the more accurate trun- 
cated expansions, is given in terms of the magnitude of a 
suitably scaled dimensionless temperature. The scaling 
factor, which is introduced in section 3, takes the 
material parameters and sphere size into account. When 
these parameters have values for which the truncated ex- 
pansions are not accurate, then the integral expressions 
for the thermal parameters can be used as a basis for an 
iterative solution of the calorimetric problem. These in- 
tegral expressions could also be used to analyze a 
calorimetry experiment conducted iu the transient 
regime, assuming that the knowledge of the surface 
temperature of the sample includes an interval (on the 
order of the characteristic decay time) that precedes the 
commencement of the cooling observations. For exam- 
ple, if a sample were held in a constant-temperature oven 
(of known temperature) long enough to become isother- 
mal, and then suddenly removed to commence cooling 
which was observed for a time interval on the order of 
the characteristic thermal decay time, the integral ex- 
pressions for the thermal parameters could be used to 
analyze the data. 

Although the analysis of this paper was carried out 
with radiative cooling in mind, only in the solution of the 
predictive problem in section 3 is the radiative cooling 
law invoked. In the analysis of the calorimetric problem, 
the cooling law is never specified. AH that is assumed is 



that the time dependence of the time-rate of total energy 
loss by the sample is known (as is the time dependence of 
the surface temperature). 

The most obvious limitation of the analysis of this 
paper (aside from its restriction to spherical symmetry) is 
the exclusion of the possibility of phase change. That is, 
the spherical sample is assumed to be either entirely solid 
or entirely liquid throughout the experiment. In addition 
to this limitation, the analysis incorporates two approx- 
imations, the more significant being the neglect of the 
spatial variation of the diffusivity in the interior of the 
spherical sample. That is, the diffusivity is assumed to be 
a function of the surface temperature (which is a func- 
tion only of time) rather than a function of the interior 
temperature (which is a function of the radial coordinate 
as well as of time). It is shown in section 2 that this ap- 
proximation amounts to neglecting a very small term in 
the Fourier equation that has the form of an effective 
heat source, but, as explained in section 5, this effective 
heat source can be taken into account (if necessary) by a 
simple iterative procedure. The other approximation, 
whose effect is completely negligible, is the neglect of the 
spatial variation of the mass density of the sample. That 
is, the overall change in average density with 
temperature is taken into account, but at each instant 
the density throughout the sample is assumed to be 
spatially constant. In other words, as in the case of dif- 
fusivity, the density throughout the sample is assumed to 
be a function of the surface temperature rather than of 
the interior temperature. 

The literature relevant to predictive solutions of the 
Fourier equation is old, vast, and still growing [3,4], 
However, this literature is almost exclusively devoted to 
the initial-value approach to the problem which requires 
that at some instant the interior temperature field must 
have some exactly specified form (most commonly, a 
given uniform temperature). This point of view, 
however, is physically inappropriate to the calorimetric 
problem because what is usually known is the history of 
the environment to which the sample has been exposed 
(i.e., the history of its surface temperature), and not the 
interior temperature field at any instant. It is true that, if 
the sample is kept in a constant-temperature oven long 
enough, its interior temperature will indeed be spatially 
uniform, but this is a special case. It would be physically 
more natural to replace the initial specification of the in- 
terior temperature field with the specification of the sur- 
face temperature history back to t = — °°. (It is shown in 
section 2 that as a practical matter it is only necessary to 
know the surface temperature during a very short period 
of the past.) There is a well-known solution to the 
Fourier equation (cf. for example p. 247 of Ref. [3]), 
which has the form of a convolution of the surface 
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temperature with the well-known diffusion kernel, but 
this solution is inappropriate to the calorimetry problem 
because it has a singularity at the center, and so (in the 
absence of a point heat source) can only be used to 
deseribe the temperature field in an infinite medium sur- 
rounding a spherical cavity. 

The mathematical literature that is directly relevant to 
the determination of the thermal parameters from the 
observed surface temperature [5-13] unfortunately has 
remained within the framework of the initial-value ap- 
proach. Because an arbitrarily specified time- 
dependence for the surface temperature is generally in- 
consistent with a previously specified initial interior 
temperature field, the problem is over-specified, and cer- 
tain compatibility conditions must be satisfied before the 
problem is well-posed. The derivation of these conditions 
has been an important theme in this literature. (The 
whole question of compatibility becomes irrelevant, of 
course, when the surface- driven solution is used as the 
basis of the analysis.) The thermal parameters have been 
expressed most commonly as the solution of an integral 
equation, but the most general case considered so far has 
allowed only one of the two parameters to be an 
unknown function of temperature, the other being an 



unknown constant. Because these solutions are very dif- 
ferent in form from the expressions given in this paper, 
and because the geometry considered was planar (either 
slab or semi-infinite medium) rather than spherical, no 
attempt has been made to compare the results of this 
paper with the previously derived expressions for the 
thermal parameters. 

1.1 Notation 

The main analysis involves dimensionless quantities, 
which are designated by bare letters, whereas the cor- 
responding dimensional quantities are indicated by an 
asterisk. Time-independent unit quantities (also dimen- 
sional) are indicated by a caret. The relations existing 
among the three types of quantities are given in table I 
which also serves to define most of the notation. {A few 
more symbols will be introduced as needed.) Table I also 
shows how the various dimensional quantities depend on 
the radius of the sphere. Because of thermal expansion, 
both the dimensional radius R* and the dimensionless 
radius R are variable, but the unit radius R is an ar- 
bitrarily chosen constant. The R -dependence of the 
various dimensionless quantities has been defined in 



TABLE I. Basic Notation for Calorimetric Problem . 



Dimensionless 


Dimensional 






Quantity Symbol 


Symbol 




Remarks 


Radius of sphere R it) 


R* = RR 






Radial distance r 


r* = rRR 




0<r< 1 


Gradient operator 7 


V* = RRV 




• 


Mass density <j(tl 


q* = qq/R* 






Q =constant 


Specific heat c(0) 


c* = cc 






' k* = Q*a*c* 


Thermal conductivity A: (9) 


k* = (k/R)k 






k = Qac 


Thermal diffusivity a{&) 


a* = R*aa 






k = Qac 


Total hemispherical 








emissivity t{T\ 


e* = el 






Linear time t 


t* = ti 




t = R 2 /a 


Nonlinear time tit) 


dt = aiTHt = 


a*dt*/R* 2 




Surface temperature T{t 1 


T*= TT 






Interior temperature 0(f,r| 


0* = 07" 




(0) r=1 = T 


Surface specific enthalpy hit) 


h* = hh 




h =cT 


Interior specific 








enthalpy n(f,r) 


i* = n^ 




. K=i = * 


Total mass 


M*=|"R*V 


= (|"» 3 e)c = 


= Mq 


Total enthalpy H 


H* = M^J = 


= \McT)qH 


H = n av 


(average specific 








enthalpy) 








Independent Reference Quantities: R, q, T, c, a 


(all constant) 







515 



such a way that the dimension less Fourier equation is 
completely independent of the effects of thermal expan- 
sion. (Cf. Sec. 2.) During an observational run, ./?*((*) 
can be measured optically along with T*(t*), and then 
R(T*( can be calculated. When, for example, this RiT*) 
is entered into the expression given in table I for a*, the 
contribution to the T*-dependence of a* that results 
from thermal expansion is automatically taken into ac- 
count. Throughout the paper, except for section 3 that 
deals with the predictive problem, it will be assumed that 
R*{t*), T*U*), and dH*/dt* are given functions of t* 
resulting from the experimental observations. From 
these the dimensionless functions Rft), Tft), and H[t) = 
dH/dt can be directly calculated, so it will be assumed 
that these too are given functions. An overhead dot will 
indicate differentiation with respect to the dimensionless 
linear, time t. It should usually be possible to choose the 
unit time f to be a convenient multiple of some ex- 
perimentally defined time interval, such as the interval 
between ohservational readings. Differentiation with 
respect to the nonlinear time coordinate t will be 
designated as follows: dk/dr = H^K Although it will 
often be desirable to choose q, a, and c to be close to the 
values of q*, a*, and c* at T* = T (which means that 
the corresponding dimensionless quantities will be close 
to unity at the reference temperature T), this is not 
necessary. 

2. Surface-Driven Solution 

If the spherical sample is imagined to be immersed in 
a heat reservoir of variable temperature, then changes in 
the interior temperature field are driven by the pre- 
scribed changes in the surface temperature. Assuming 
the absence of any interior heat sources, it follows that 
the interior temperature field is uniquely determined by 
the past history of the surface temperature up to the 
present moment. In mathematical terms, this cor- 
responds to the "particular" or "driven" solution of the 
Fourier equation, with the surface temperature playing 
the role of the "driving functiorj." This is not the most 
general solution, because it does not include the 
homogeneous solution which describes the decay of an 
arbitrarily specified initial interior temperature field. It 
is well known that the most slowly decaying term in the 
homogeneous solution has a time dependence propor- 
tional to exp(— n 2 r) where t is the dimensionless time 
measured in the natural time unit defined in table I. 
Neglecting the homogeneous solution amounts to assum- 
ing that the interior temperature field has been subjected 
to no influences other than its external environment for a 
period of time t that is long enough so that 
exp{— 7i 2 ()«l. 



The time-rate at which the sample exchanges energy 
with its surroundings is determined by the history of the 
surface temperature up to the present moment. In fact, it 
is just equal to the time derivative of the total interior en- 
thalpy of the sample. Thus, once the time history of the 
surface temperature has been specified, the time-rate of 
energy loss or gain of the sample is completely deter- 
mined. The analysis of this section leads to expressions 
(summarized in tables V & VI) relating the time-rate of 
total energy change of the sample to the surface 
temperature (or more exactly, the specific enthalpy at the 
surface), and these expressions suffice for the analysis of 
both the predictive and the calorimetric problems. 

The dimensional Fourier equation is given in the two 
forms (la) and (lb) of table II, the only difference being 
the representation of the part of the heat flux that results 
from radial motion caused by thermal expansion or con- 
traction. In eq (la) it is represented in terms of the 
material velocity v* at a point r* that is fixed in the 
laboratory (inertiall frame, whereas in eq (lb) the motion 
is taken into account by the fact that the time derivative 
is taken with respect to fixed r rather than fixed r*, 
where r is the dimensionless radial vector that is attached 
to a particular material particle and moves with it. 
Although v*, which is the material velocity associated 
with thermal expansion or contraction, is negligibly 
small, the point to be made is that the right-hand side of 
the dimensionless Fourier equation given in eq (2) is 
rigorously correct, and the fact that the time derivative is 
taken at constant r rather than r* does not represent an 
approximation. 

In eq (3) the internal enthalpy density 17 is introduced 
in order to replace the internal temperature 0. This 
replacement is doubly advantageous: First, a com- 
parison of eqs (2) and (4) shows that it reduces the 
number of thermal parameters that appear in the equa- 
tion. Second, the enthalpy density is really the quantity 
of physieal interest, because the objective of the analysis 
is to integrate it over the volume of the sample in order to 
arrive at an expression for the time-rate of change of the 
total enthalpy (heat content) of the sample. 

Equation (4a) still contains the temperature- 
dependent diffusivity a(0), and this fact not only com- 
plicates the equation, but also prevents it from being 
universal in the sense of having the same form regardless 
of the material properties of the sample. If the diffusivity 
were a function only of t and not of r, it could be 
eliminated from the equation by replacing the linear 
dimensionless time t with the dimensionless nonlinear 
time t as indicated in eq (5b). In fact, this device for 
eliminating the diffusivity has been used before [5,6], 
The same device would also eliminate a from the equa- 
tion if it were a function of the surface temperature T(t) 
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TABLE IL Basic Differential Equation . 



Eg. No. 



m 

(4a,b,c)r 

{5aM 
|6a,ht 

m 

iSaM 
I8c| 



V *•(&* V *&*} = e *c*[3 Q*/ 3t*) E » iv*'V *&*]. = q*c*{ 30*/ 3t*r r 
\T«(/c¥0)=ec(36/34 

<f ^ = c{qu e 

¥■(!<» VijJ — (3^/3r)^ a — k/fQc; q = constant 

V *3 + 7 • ir = On/ 3t), ; dx = a(Tfe])<fc = a\t)dt 

uS ^ t _ am lVn; v. u = qell 

3 2 |rn} 3frn) ^ _ _ I £_ f p, °< e > , , a 1 
Sr 2 3t 



-^ = r l f r {[i-W]' 2 -a7l=° 






rather than of the interior temperature Q('r,f )'. Even in 
the latter case, however, introduction of a nonlinear time 
t baaed on afTfi])' succeeds in eliminating most of the 
a-dependence from the equation , as shown hy eqs (5a). 
and (<6a)'. What remains in the equation Is a very small ef- 
fective heat source density term q e£f which has the form 
of the divergence of an effective heat flux u defined in eq 
(6a). This heat flux vanishes at the surface of the sampler 
and for this reason, when Gauss' theorem Is applied to eq 
(Sal In order to arrive at the expression given In eq (8h( 
for the- time-rate of total enthalpy change, the term In- 
volving u makes no contribution. Therefore, because eq 
(!8W leads to the equation (eq ( 1 91 of table V) from which 
the rest of the analysis follows, it Is evident that, at least 
to first order, the introduction of the nonlinear time T has. 
succeeded In reducing, the problem to the solution of the 
universal equation that results if q e f f = & In eq (7)L The 
analysis of this paper Is based on this approximation- If 
more accuracy should be required, then the solution for 
ijfjyrF that Is given In eq (9) or eq (12). of tahle III could 
be substituted into the right-hand side of eq (7), and an 
additive correction to n; could be found which In turn 
would lead to an additive correction to H which could be 
Introduced Into the calorimetrlc equations of section 4. 
The way this would be done is explained in section 5. 

The solution to eq (IK (with on the right-hand side) 
that satisfies the boundary condition, stated in. eq t8a)i can. 



be written In the form of eq (9} In table III- The 
polynomials p n (rt are characterized by the property 
stated In eq (10a)), and can be generated by successive In- 
tegration. The first four polynomials are given In eq (II), 
and are plotted In figure 1 . The fact that eq (9) does in- 
deed satisfy eq (7)' (with 0: on the right-hand side) can he 
directly confirmed using the property stated in eq (10a). 

As indicated In eq (12), the solution can' also he ex- 
pressed In terms of the odd-order Bernoulli polynomials 
H^-tl^ where x = M(l-r). The properties of these 
polynomials that are necessary to verify that eq (12)! Is in- 
deed the desired solution of eq (7)' are stated In eqs (13) 
and (14). (See for example,, pp. 804-811 of Ref. [14] or 
pp. 19, 25-29 of Ref. [I5].jt 

An explicit expression for H"' = dH/dr In terms of 
£.(n) = dph/gft can ne derived by substituting; eq (1 2) Into 
the left-hand side of eq (8b) and using the relation stated 
in eq (16)1 of table IT between the even-order Bernoulli 
numbers B^ and the Rlemann Zeta function £(2n)'. The 
resulting relation is given in eq (1 9J. of table V. Using the 
numerical values for £(2n)i that are given in table IV, eq 
(20) results, which can be then Inverted to yield eq (21), 
which will play an important role In section 4. 

The expansion given In eq ( .19} assumes that h(rt is an 
analytic function all of whose derivatives exist. If In ad- 
dition it remains finite for all t, it can he shown that eq 
tl9t> is, equivalent to the integral equation given In eqs 
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Table 1 1 1 . Surface-Driven Solution for Interior Enthalpy Density . 



Eq.No 



CO 

(9a,b) n(r, T ) = h{r) + I (-1)" p„(r) A<"> (t) ; h M = d n h/di n 

n = l 

<10a,b) <f 2 (rp„)/<fr 2 = - rp^ or V 2 p„ = - P„-i 

(11a) Pl = (l-r 2 )/6 

(lib) p 2 = (7 - 10r 2 + 3r*)/360 

(lie) p 3 = (31 - 49r 2 + 21r 4 - 3r 6 |/15,120 

(lid) p 4 = (381 - 620r 2 + 294r 4 - 60r 6 + 5rVl, 81 4,400 



► Cf . Fig. 1 



(12a,b) n (r, T ) = M T )-il 

' « = 1 



22n + l 
i.i - 



Ba,+ iWfc W W; *sl(l- r > 



(13a,b) [Ba,+iWUo = ^+1^ = ° = P&+1WU1 = *j»+i«>> = 

(14a,b) dB ^M = _i (jjn+ijB^w. d 2 *2n+I<*> = ^JJn+lH&ilBj^W 



<fr 



dr 2 



(22a) and (23) of table VI. The kernel V of the convolu- 
tion integral defined in eq (23) is an effective memory 
function that weights the very recent past most heavily 
and totally forgets events that happened more than half 
a natural time unit in the past. This memory function is 
defined by eq (24) and is plotted in figure 2. Its argu- 
ment is defined by eq (27a), and as shown in eqs (271> 
and c) can be expressed in terms of a(t) and the dif- 
ference {t~t' ) between the present time t and some past 
time t' . Figure 2 shows that for £ less than 0.1 natural 
time units the simple function f" defined by eq (25b) is 
essentially indistinguishable from I", For larger £, the 
first term in the summation of eq (24) should serve to 
represent l~ with sufficient accuracy for most purposes. 
As eq (26) indicates, the normalization of V is such that if 
hS^ is constant, then the retrospective weighted average 
h^ defined by eq (23) will Just be equal to h^\ If, 
however, hS^ varies drastically during half a natural time 
unit, which could be the case when a sample first starts 
to cool, then the weighted average /i l2l U) will differ 
markedly from the instantaneous value /t ,3l ((). In such a 
case the integral eq (22a) will be more aecurate than the 
equivalent truncated expansion given in eq (20). The 
series expansion (19) can be derived from the integral 
equation defined by eqs (22a) and (23b) by expressing 
A' 2 '(t — £) as a Taylor expansion about t, integrating by 
parts, and making use of the definition of £(2n) given in 



eq (15) in table IV. 

Finally, it sho uld b e noted that it is evident from eq 
(23b) that when /i (2 '(t) is differentiated with respect to t, 
the differentiation can be taken inside the integration, 
from which it follows that eq (22b) results from differen- 
tiation of eq (22a). Obviously, an infinity of such equa- 
tions can be generated by repeated differentiation. 



TABLE IV. Riemann Zeta Function for Even-Integer Argument . 



Eq.No. 



(15) <(2n) s I m" 2 " 

m = l 



(16) 



*& s [Ba,Wl=o = 



(-1Y>-H2n)\ c,(2n) 



2 2n-l 



„2n 





2n 


daii/n 2 " 


<S2n) 


2 


1/6 


1.64493 


4 


1/90 


1.08232 


6 


1/945 


1.01734 


8 


1/9450 


1.00407 


10 


1/93,555 


1.00099 
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TABLE V. Total Rate of Energy Loss: Differential Representation . 



Eq. No. 



(17a,b,c) /f (1 » = dH/dr ; H = dH/dt ; H™ = H/a\t\ 

(18a,b) A (1 » = cMrUH-r) = c(t\t{tV a(t\ ; ft (n > = d"h/dT» 

(19) /f n >(-r) = - 6 I <-D n k(2n)/n 2n ] A<»>(t) 

n = l 

{20 ) mi = A* 1 * - is * m + 313 * (3) - jk h ™ + ■ ■ ■ 

(21) A m (T) = fl <l> + Xtf< 2 >-^Hl3> + _^H<4, + . 



Eq. No. 



TABLE VI. Total Rate of Energy Loss: Integral Representation. 



<22a,b) H™ = h™ -^ A (2) . #12) = ft (2( _ _L A (3) . etc . 
(23a,b> A< 2 » s J T |-(t-t' ) A< 2 >(t' ]dz' = J™ H£) h™{r-t.\dl, 



(24) r(S) a 90 2 M" 2 expHmnR] 

m = l 

(25a) T(0 = T(Q for < £ < 0.1 
<25b) HO e 15[1-3(2n/^7S -0] 
(26a,b> J T Ht-t' Wr' = J V{QdZ = 1 

-co 

(27a) £ h t-t' = S\ a\t')dt' 



(Cf. Fig. 2) 



t' 

00 



(27b) £ia; t-t') = l 



(-D- 



icPa/ dt*)it~t • ) n 



n=0 (n+1)! 
(27c) £<<*; f-f ' ) = ait)it-t' ) - \ <Ht)(t-f ) 2 + \ dit)it-f ) 3 - . . . 



3. The Predictive Problem 

For the purpose of estimating radiative cooling times 
and the relative magnitudes of the terms in eq (21), from 
which the calorimetric equations of section 4 are derived, 
it is useful to solve the post-transient predictive problem 
for the case of constant parameters. In such a case the 
simplifications indicated in eqs (28a-d) of table VII oc- 
cur. All of the equations of tables V and VI are still valid. 
In particular, H = if' 1 ' must satisfy eq (19), but in addi- 



tion it must satisfy the Stefan-Boltzmann radiation law 
which means that the left-hand side of eq (19) must be 
replaced by T 4 multiplied by a proportionality constant 
involving the Stefan-Boltzmann constant o* = 5. It is 
easy to show that if the unit temperature T is defined as 
shown in eq (29b), the proportionality constant on the 
left side of the specialized form of eq (19) will by unity, 
with the result that the equation has the form given in eq 
(31). A significant feature of this equation is that it is 
universal in the sense that it applies to spherical samples 
of all sizes made of any material whose thermal 
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-p,* 10 ' 



-P 3 x10 J 



P 2 x10' 



Pi x1 ° 




r(o)=i5 




FiGUBEl. DimensionleBs polynomials p (rl defined in eq8<]la-d) 



FIGURE 2. Dimensionless kernel r(£l "in convolution integral delined in 
eq (23), and approximate kerna] f {£] defined in eq I25b), 
both as functions of the dimensionless argument £= x^ t'.. 

TABLE V1L Equations for Predictive Problem. 



Eq. No. 



(28a,b,c,d) c=a=k = R=iL=Q=\* t = f; 13 = 0; h = T 



(29a,b) 
(30a,b) 

(31) 
(32) 

(33a,b) 

(34a,b) 

(35a,b) 

(36a,b) 

(37) 

(38) 

(39) 



T= T*/T; T= (fc/i?o£) 1/3 = (fc*/R*o*£*) iy3 = ( e *a*c*/J?*<j*.E*) 1/3 
H = dH/dt = -3T i ; (3rj/9r) r=] = (d6/dr) r=1 =-T 4 ' 



T 4 = 2 I <-l) n [^(2n|/n 2n ]d n r/dt" 



3 A 45 dr 2 



1 d A T 



+ 



i d 4 r 



945 dt 3 4725 dr 4 



dr/dt = -3r 4 [i- |r 3 +^r° + ^r 9 +...] ; d#/<fr = -jr* 



d 2 T/dt 2 = 36r 7 [i-^r 3 + f r 6 + ...] 



25 675 

25- 



d*T/dt s = - 756 T™ [1 - if 7* + . . . ] 

d*T/dt 4 = 22,680 T 13 - . . . 

tfTjl = | (7^ 3 -1) — ±ln r + f 5 (i - r 3 ) 



d 2 ff/& 2 = 36T 7 [l- 1 3* + ^r 6 + ,„] 

d*H/dt 3 = -7567™ [MjSj T 3 + ...J 
d^H/dt* = 22,680 r 13 - . , . 



to, = pi + 9f ) _ 4_ in (i + 9t) - i| 1T ^ 1 . r" 3 
rw = [- i t (r) + ^ J" r<fi r((-4id£] 1/4 



(Cf. Fig. 3!) 
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parameters are constants. It is possible to solve this 
equation by truncating it at the fourth derivative as 
shown in eq (32), and solving the truneated equation for 
dT/dt in terms of T by starting with the approximation 
dT/dt = — 3T 4 and iterating until a self-consistent set of 
expressions for the first four derivatives results. These 
expressions are given in eqs (33a-36a). By differentiating 
eq (30a) and making use of eqs (33a-36a|, the expres- 
sions for d n H/dt 11 given in eqs (34b-36b) can be derived. 
The 8 expressions in eqs (33-36) will be used in section 4 
to estimate the range of validity of the calorimetric for- 
mulas derived there. Numerical estimates indicate that 
these expressions are accurate to within 1 % so long as T 
<14. For larger 7\ it would be necessary to include 
higher order terms in the expansion given in eq (31 ), and 
the numerical estimates indicate that for T > 3/4 the 
convergence is so slow that this expansion has no prac- 
tical utility. Correspondingly, the calorimetrie equations 
derived in seetion 4 that are based on eq (21), which is 
derived from eq ( 1 9), cannot be expected to be accurate, 
even in a post-transient experiment, if the dimensionless 
surface temperature T based on the unit temperature 
defined in eq (29b) is larger than Vi. If the sphere size 
and thermal parameters are sucb that T > Vz, then it will 
be necessary to use equations based on the integral equa- 
tion defined by eqs (22a) and (23). In order to give a feel- 
ing for what sphere sizes and which materials will satisfy 
the condition T < ]/i, the dimensionless temperatures 
■^melting corresponding to the respective melting poults of 
tungsten (3650 K) and uranium dioxide (3150 K) are 
given in table VIII for sphere radii that approximate the 
upper and lower limits that would most probably be con- 
sidered for radiative-cooling calorimetry experiments. It 
is evident from this chart that for most praetical post- 
transient experiments, it should be possible to use 
ealorimetrie relations derived from tbe truncated expan- 
sion given in eq (2 1 1. Only in the ease of a large sample 
(R ~ 1 cm) of a poor thermal conductor (such as uranium 
dioxide) might it be necessary to use an integral relation 
in order to analyze the results of a post-transient experi- 



T ABLE VIII. Representative Values for Unit Time and 
Unit Temperature. 





k 


Tungsten 


Uranium Dioxide 




0.1cm 


1 cm 


0.1cm 


1cm 


melting 
J melting 


0.03 sec 
37,000 K 
3650 K 
0.1 


3 nee 
17,000 K 
3650 K 
0.2 


2 sec 
8100 K 
3150 K 
0.4 


200 bcc 
3800 K 
3150K 
0.8 


u 


it 2 /* 


f=i*/*a£) 1/3 = (5ac/.R6£| 1/3 





O) 



— 0.1 T 




logt 

FiGUHE 3. Universal temperature-time curve for post-transient 
radiative cooling. 

ment. It ought to be noted, however, that this might be 
neeessary even in the case of a smaller sample of a better 
conductor if the observations are based on a transient- 
type experiment. 

A universal post-transient cooling curve ean be de- 
rived by integrating eq (33a). The result is the expression 
for t(T) given in eq (37). This ean be inverted to yield the 
expression for Tit) given in eq (38). The cooling curve 
corresponding to these expressions is plotted in figure 3. 
It is evident from this curve that the slope of In T versus 
In t is almost, but not quite, constant. In fact, this slight 
variation in slope is related to the thermal condnetivity 
of the sample. It can be shown that for initial and final 
temperatures T-* and 7V* 



(40 a) 




(40 b,c) 



This equation cannot be used for determining thermal 
condnetivity from observation of post-transient radia- 
tive cooling, because it assumes that the specific heat is 
constant throughout the cooling, whereas in all pro- 
bability the In T* versus In t* curve for a real sample 
would have much more curvature than the one shown 
in figure 3, and most of this curvature would be caused 
by the temperature dependence of the speeific heat. The 
real significance of eq (40) is that it (together with Fig. 
3) illustrates how difficult it is to make a reliable deter- 
mination of thermal conductivity (or diffnsivity) from 
observations of post-transient radiative cooling, 
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especially when these observations are made for 
temperatures near the lower end of the curve shown in 
figure 3. In contrast, it is very easy to determine the 
specific heat in this range since to a good approxima- 
tion it will be given by c =H/T. For an accurate deter- 
mination of thermal conductivity it will probably be 
necessary to use a sphere that is large enough so that 
the dimensionless temperatures involved fall well above 
those shown in figure 3. In such a case it would be 
necessary to analyze the data using the integral expres- 
sions given in section 4, rather than the truncated ex- 
pansions. 

It should be noted that when the expressions given in 
eqs (33a-36a) are substituted into eq (9a) taking the 
simplifications stated in eqs (28a-d) into account, a com- 
plete solution for the interior surface-driven solution in 
terms of the surface temperature T results. If the expres- 
sion for T(t) given in eq (38) is substituted into this, an 
explicit expression for the interior temperature field 
0(r,r) results. If it were desired to extend the validity of 
this solution to values of T larger than Vi , this could be 
done by using the integral equation given in eq (39) in 
table VII as the basis for an iterative solution. Equation 
(39) was derived from the integral equation defined by 
eqs (22a) and (23), making use of eq (30a). The idea of 
reducing the problem of solving for the interior 
temperature field to the problem of solving an integral 
equation involving only the time dependence of the sur- 
face temperature is not new. It has been done for a semi- 
infinite medium with a plane surface [16]. The integral 
equation that resulted was derived from the diffusion 
convolution integral mentioned in section 1. However, 
this approach is not appropriate for the present problem 
because, as noted in section 1, the diffusion convolution 
integral represents the temperature field in an infinite 
medium surrounding a spherical cavity, rather than the 
field within a finite spherical medium. 

4. The Calorimetric Problem 

In adapting the expressions derived in section 2 to the 
problem of deducing the specific heat and the thermal 
diffusivity from observational data, the choice made for 
the unit temperature T can be arbitrary. It is not 
necessary to use the unit temperature defined by eq (29b) 
of table VII, although this choice is appropriate for the 
purposes of designing an experiment, and for determin- 
ing whether the various expressions derived- in section 
4.1 from eq (21) are accurate, or whether it is necessary 
to use the alternative integral relations discussed in sec- 
tion 4.2. 

Once it is a question of analyzing existing data, 
however, it would generally be more convenient to define 



T so that the dimensionless temperature T is close to uni- 
ty. If the thermal parameters are already known for the 
cold end of a post-transient cooling run, then it would be 
natural to choose this cold temperature as the unit 
temperature, and correspondingly the unit specific heat 
and unit diffusivity would be chosen to be equal to the 
known values of these parameters at this cold 
temperature. If the true values are not known, then 
estimates would suffice. Nowhere in the analysis is it 
assumed that these estimates are close to the true values. 
For example, if one were analyzing data for samples of 
different materials, it might be most convenient to make 
a single choice of unit quantities to be used for all of the 
different materials. 

In all of the expressions given below, Hit) — dH/dt 
and T(t) are regarded as given functions of time that 
result from independent simultaneous observations 
made by two different instruments. If a reliable cooling 
law exists and is known, then H can be expressed as a 
function of T and eliminated from the equations. In the 
case of radiative cooling, this would require that the 
temperature dependence of the total hemispherical 
emissivity £(71 be known. 

4.1 Truncated Expansions 

The calorimetric formulas, which were the principal 
objective of this analysis, are given in eqs (41) and (42) of 
table IX. The expression for c follows directly from eq 
(40), which is simply eq (21) of table V multiplied by a[t). 
The expression for a was derived from the ratio of the 
time derivative of eq (40a) to eq (40a) itself. Both expres- 
sions for c and or have the form of a power series in an ex- 
pansion parameter | = 1/15ct. The coefficients of these 
power series are functions of the four quantities defined 
in eqs (43a-d), the leading terms of which are ratios of 
different time derivatives of T and H. Equations (44) 
and (45) give approximate expressions for c and a that 
are valid in the limiting case in which only the leading 
terms in the expansions must be retained, and c and a 
are essentially constant. 

In the discussion that follows, it will be assumed that 
the unit of diffusivity & has been chosen so that for the 
data under consideration the dimensionless diffusivity a 
is of order unity. Then £ ~ 1/15. (If a different choice of 
a were made, the change in £ would be compensated by 
changes in the values of the quantities defined in eqs 
(43a-d.) Because the expansion parameter | involves a, 
and the coefficients in the expansions for c and a in- 
volve c and a, it is evident that eqs (41) and (42) must be 
solved iteratively, with the first iteration based on the 
assumption that £ = c = a = 0. The range of con- 
vergence of this procedure can be estimated by using the 
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expressions for d n T/dt n and d a H/dt n given in eqs (33- 
36) of table VII to evaluate all of the terms in eqs (41) 
and (42). When this is done the following expressions 
result: 



A =-12 7 ,3 (l-|7 ,3 + |7 ,6 + . 

B =-12 7 ,3 (l-|7' 3 + ^T 6 + 
5 2a 

C =252TMl-|^r 3 + ...) 
D = -7560T 9 + . .. 



A-B 



AD _ 



A-B 



= 12,600 T 6 + 



H/f= i+| T^-gr 6 

-i- n5 i -r- ... 



1 / C-AB 



15 I A-B 



= 1-17-3 + 



(46a) 
(46b) 
<46c) 
(46d) 

(47a) 

(47b) 
(48) 

(49) 



When these expressions and | ~ 1/15 are substituted 
into the right-hand side of eq (41a), it reduces to 1 + 
0(T 12 ) which (since the left-hand side is c = 1) is just the 
identity that is to be expected in view of the fact that eq 
(41a) is simply a reformulation of the same equation 
from which eqs (33-36) were derived. Similarly, the 
right-hand side of eq (42) reduces to 1 + 0(T 6 ). Thus (to 
the accuracy of the truncation) the leading factors in eqs 
(41a) and (42), which are now expressed by eqs (48) and 
(49), are just the reciprocals of the respective square 
brackets on the right-hand sides of eqs (41a) and (42). 
For this reason, the speed of convergence of the 
calorimetric formulas can be estimated by inspecting eqs 
(48) and {49). These indicate acceptably rapid con- 
vergence for T < Yz which, of course, is the same range 
of convergence that was noted in section 3 for the val- 
idity of the iterative solution of the predictive problem. 

In the case of the approximate limiting expressions 
given in eqs (44) and (45), all of the terms of eqs (41a) 
and (42) that involve | were thrown away, and only the 
leading terms were retained. In addition, all of the terms 
of eqs (43 a-d) involving c and a were thrown away. The 
validity of this latter approximation can be answered 
only on a case-by-case basis, but the validity of ignoring 
the terms involving | can be estimated by means of eqs 
(48) and (49) since the terms involving T represent the 
error in these formulas, because in this case the correct 



values of c and a are unity. In the chart below the 
magnitudes of the two leading error terms are listed for 
several values of T. 



4^3 

S 



0.1 



0.10 

0.013 

0.0008 



0.15 

0.019 

0.001 



This chart shows that for T < */i> the error is less than 2 
percent, and becomes increasingly smaller the smaller T 
becomes, i.e, as T enters the extreme post-transient 
regime. Even for T = V2, the error is not so large as to 
destroy the usefulness of eqs (44) and (45) for generating 
the first iterative solution for cit) and ait) which is then 
substituted into the right-hand sides of the more accurate 
formulas given in eqs (41) and (42). If a is constant, it 
follows from eq (45 d) that H cc exp[-15a Q t], so that in 
the extreme post-transient regime a Q can be estimated by 
fitting the observed function Hit) to an exponential 
decay. 

It should be noted that, in order to make these 
estimates, it has been necessary to define the dimen- 
sionless T in the manner indicated in eq (29) of table 
VII. This automatically takes the sphere size and ther- 
mal parameters of the sample into account. However, for 
an actual application of the calorimetric formulas of 
table IX, it is not necessary to do this. One may use any 
convenient scaling factor to define the dimensionless T. 
The validity of the formulas would then be indicated 
directly by the convergence behavior of the numerical 
iteration process. 

The solutions of eqs (41) and (42) are c(t) and a{t). 
However, because Tit) is known from observation, these 
solutions can be converted into c{T) and aiT), which are 
the desired expressions for the temperature dependence 
of the thermal parameters. 

4.2 Integral Relations 

If a numerical application of the calorimetric formulas 
shows that the £ 3 term in eq (42) is comparable in 
magnitude with the | 2 term, or if the | 2 term in eq (41a) 
is comparable with the | term, this is an indication that 
the neglected higher-order terms are not really negligi- 
ble, and that the calculated functions ait) and c{t) are not 
reliable. One could, of course, include higher-order terms 
in the equations, but truncation must occur at some 
point, so the net result would be only a slight extension of 
the range of validity of the equations. Moreover, the 
higher-order terms involve higher derivatives of Tit) 
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TABLE IX. Calorimetric Problem: Truncated Expansions. 



Eq. No. 



|40a,b) 



(41a,b) 

(42) 

<43a> 
(43b> 
(43c) 

(43d) 

(44) 

(45a,b) 

(45c,d) 



A=cT = J &+^|^— ^^+^4^.+ ...]; 



15 dt 



35 dt 



\hw=K h^ = -L^'' h< 3 >= 1 dHW 



dt a dt 

i3 /"i t 2 a. a n i «3 



] 575 dt 
] 



c = (H/T)[l + S^-(f C)^ 2 + (fl»£ 3 +...]; |= l/15a 



= L,C-^Bhi _ 3,£»-^C, f2 _2 , ^£> i ,3 , 



.<* 



^> [1 - ^ 



-) I 2 -= <- 



] 



^ s d In MD/d, = d In {cf/a)/dt = (f/t) + [(c /c ) - (a/a)] 

fl = d In H^/dr = d In (tf/a)/dr = [H/H) - {a /a) 

C = a^H^/H^ = {H/H) - 3(a/a) + [(a/a) + Ca/a) 2 ] 

D = »3fl1«/jj1l) = (#/#( - 6(a/a) (tf/if) + [15(a/a) 2 - Ma/a)]{H/H) -[fa/a) + 15(a/a) 3 

- lOtaa/a 2 )] 



c =(c) |=( ; =i=0 =/f/r 

1 ,C - AB X 



i [TH-ra 



_ / i i ii-t s*£>\ i \i n — i n 



1 

15 



r i? 






15 dt j< dt t 



Table X. Calorimetric Problem: Integral Relations. 



Eq, No. 



(50) 
(51) 
{52) 

(53) 

(54) 
(55) 
(56) 
(57a,b) 



ait) = { [(h-H)/(h -H)] - (a/a) } [A< 2 »/ft< 3 >] 

ft (r) = if + \ t { [(h -H)/(h - H)\ - (a/ a) } [{$%)* /W\ 

hi3i=r ri-f i J^^ld,' 
-'-m d( |^ <r(r ) rfi o(r ) J J 

r = 90 2 (mn)- 2 exp[- (mn) 2 £] 

m=l 

£ = /' a{t")dt" = a(t)(t-f ) - ^tHt-t 1 ) 2 + faWt-t 1 ) 3 - , 

c[t)=h{t)/f(t) 

c(T> = citffl) ; a(T) = a(f[71) 
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and HitU and the error involved in .extracting these from 
the experimental data becomes ever greater the higher 
the order of differentiation. For these reasons, it is better 
to use an Iterative procedure based -on eqs (5fl( and (51) 
of table X which involve the integrals defined in eqs (52) 
and (53 L These equations were derived from eqs (22a) 
and (22b) of table VI. Because h(t) = cit) 7%), eq (51:) is 
really an equation lor.cftk but it is simpler to regard hit) 
as the unknown function and, after this has been found, 
to invoke eq 456) to find c(r), and eqs 157.a,b : ) to find c{Ti 
and aVD. The kernel F of the integrals is defined by eqs 
(54) and (55), but the approximation based on eq 125b) 
and discussed at the end of section 2 would simplify the 
calculations. Inasmuch as the temperature dependence 
(and hence the time dependence) of a is usually weak, 
and las Fig.. 2 indicates) T will usually vanish in a time 
interval that is short compared with the time required for 
a to change by a significant amount, in all but the most 
extreme of transient experiments it would be justified to 
drop all but the first term in the expansion lor £ given in 
eq (55). 

The iteration could be started with an -Mt) calculated 
from eq (45) and an Mi) = c(i) T\t) where c)t) is found 
from eq (44). These approximate functions would be 
substituted into the right-hand sides of eqs (50) and (51), 
which would yield new (presumably improved) approx- 
imations. QnestionB of convergence and numerical 
stability of this procedure have not yet been investigated. 

5. Discussion 

The foregoing analysis took the temperature 
dependence of the specific heat fully into account, but 
the interior spatial variation of the diffuslvlty was 
neglected. This amounted to neglecting an effective heat 
source density in the Fourier equation, but it was 
pointed out in section 2 that this neglected term could be 
taken into account in an iterative fashion. This would 
give rise to an additive correction d a r\ to the interior en- 
thalpy density field. It was noted in section 2 that when 
Gauss' theorem is applied to the Fourier equation, the ef- 
fective heat term makes no direct contribution to the 
resulting equation (eq (8b > in table II). It does make an 
indirect contribution, however, in the sense that it pro- 
duces an additive correction to the radial derivative of 
the interior enthalpy field. Thus, eq (8b) must be re- 
placed by 

[ 3 .(r, + d ij1 ,/ 8 r]^ 1 = ^}(H + 6JQldT (58a) 

where H is given by eq (8c) and (J^ff is given by 

<Lr#= 3 X i6 a n\r 2 dr . (58b) 



The total time-rate of energy loss of the sample, which is 
to be identified with tbe observed energy flux, is given by 
H tot = H + 6 a H. It is important to note that, in making 
the correction to the calorimetric formulas of section 
4, fl toJ must not be substituted in place of H. The reason 
for this is that these formulas were all derived from eq 
( 19) (or its integral equivalent given in eqs (22a) and (231), 
which in turn was derived from eq (12)1, which is a rela- 
tion between the uncorrected Interior enthalpy field and 
the time-dependence of the surface enthalpy, which is 
nnaff ected by the correction because the boundary con- 
dition stated in eq tSd) continues to be valid. Thus the H 
that appears in all of the formulas of section 4 must con- 
tinue to refer to the average value of the uncorrected in- 
terior enthalpy field, which means that the right way to 
make the desired correction is to substitute the right- 
hand side of 



H — Hiot^iJ 1 



(59) 



wherever H appears In a formula, and to Identify H tot 
with the observed heat flux. 

The calorimetric formulas of section 4 yield c[t) and 
ait) as continuous functions of time from which c(T) and 
a(T) are found. If, however, spline representations of 
c{T) and a(T) are used, then the unknowns are the spline 
coefficients, which are constant numbers. Expressions 
for the coefficients as weighted Integrals involving the 
observed functions T{i) and Hit) could be derived from 
either the expansion or the integral forms of Xhe 
calorimetric formulas. Because the spline coefficients are 
expressed as integrals of the observed data, there would 
be an automatic smoothing, which could be ad- 
vantageous in the case of noisy data. 

Although the calorimetric formulas derived In section 
4 were Intended to be used with observational data from 
a single observational run using a single sample, It would 
also be possible to use them with data from two different 
runs over the same temperature range using a large and 
a small sphere of the same material. The radius of the 
small sphere would be made small enough so that the 
observations would be in the extreme post transient 
regime (T < l /i) where the accuracy of eq (44) for c 
would be good. The functional dependence for c{T) 
found in this way conld then be substituted into the left- 
hand side of eq (41a), and data from the run with the 
larger sphere could be inserted into this equation, which 
would be solved for a. This two-sphere approach had 
been suggested by Colwell [1,2] when he first proposed 
radiometric calorimetry of freely cooling spheres. 

Because the calorimetric formulas have been derived 
from an analysis that did not require a knowledge of the 
cooling law, but rather only the time-dependence of the 
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total rate of heat loss (or gain), they could be applied to 
any situation in which the interior temperature profile of 
a sphere is determined by the changing temperature of 
its external environment. For example, in the case of dif- 
ferential scanning calorimetry, using the calorimetric 
formulas would permit a determination of the thermal 
parameters of a spherical sample even when the time- 
rate of change of the surface temperature of the sample 
(i.e., the scanning rate) was so fast that the sample in- 
terior would be far from isothermal. This would permit a 
faster scanning rate, which would in turn cause larger 
heat fluxes which could be measured with greater preci- 
sion than the small ones that result when the scanning 
rate is slow enough to keep the sample interior essen- 
tially isothermal. Moreover, the differential scanning 
technique would no longer be limited to the measure- 
ment of heat capacitance, but could also be used for 
measurements of thermal diffusivity. 

Finally, because the analysis does not assume that the 
measured quantities are monotonic in time, it could be 
adapted to modulation calorimetry in which the sample 
surface is subjected to a periodically varying 
temperature and the magnitude and phase lag of the 
heat flux as a function of the frequency of the 
temperature variation are the measured quantities from 
which the thermal properties are deduced. Although the 
basic approach of this paper would still be applicable, it 
would be necessary to subject the surface-driven solution 
of section 2 to a Fourier analysis in order to express the 
various quantities as functions of frequency rather than 
of time. 



It was noted in the introduction that the analysis of 
this paper was the outgrowth of a proposal first made by 
Dr. J.H. Colwell of the National Bureau of Standards, 
and throughout the course of this work the author has 
been the beneficiary of frequent very helpful conversa- 
tions with Dr. Colwell. 

This work was supported by the Materials Processing 
in Space Division of the National Aeronautics and Space 
Administration. 
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